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ABSTRACT 



The advent of tilt rotor technology asks for rotors that 
have different twist and RPM requirements in hover and in 
forward flight to optimize for operational conditions. In 
order to get an assessment of the capabilities to fulfill 
these requirements this report presents a mapping of twist 
angle variation as a function of RPM and laminate orientation. 
The basic laminate for the six models as well the D-shape spar 
that represents the structurally active part of the blade is 
assumed to be constant ( 0®/90^/0®/90®/90®/0‘^ ) . This six layer 
cross ply laminate is chosen as it provides the necessary 
extension-twist coupling without a hygrothermally induced 
twist that is highly undesirable. The couplings and trends in 
the models are visualized in carpet plots, one for each model, 
in an attempt to establish a method to answer the basic 
question of the magnitude of twist angle available due to a 
particular geometry, material and load system. 
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A. OVERVIEW 



In the design process of composite material structures it 
is often noticed that the best characteristics of the material 
are not used. The normal procedure is to tailor the material 
to approach the properties of isotropic materials, avoiding 
the ’’undesirable” couplings between extension and twist foi‘ 
example . 

That traditional approacli worked for lielicoptei design. 
However, iii the tilt rotor, the idea of changing the twist 
angle of' tlie blade to oi>timize for either hover or forward 
flight came into the scenario. 

This study will use a blade spar to study the sensili\ity 
ol' the ext ens i on- tw i s t coupling due to changes in RPM , la>er 
angles definitioir, changes in laminate orientation and also 
the effects of each of these configurations in the first four 
natural \ibration modes of the blade. 

In the analysis, a program named CASA/GTFTS (Computer 
Aided Structural Analysis/Graphical Interactive Finite Element 
Total System; is adopted. The program makes it easier to 
xisualiz^' tlie effects of the coupling in the structure, antj is 
^ecn to an eff icieiit design tool. Tht i*esuJts obt.Ained 
through the filiate t-Jemont pr-ccedurc' i'oi eacii ianuntte 
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configuration, in a total of six models will be presented in 
3-D graphs that show the trends for each model for the layer 
angle (top layer as a reference), the RPM and the twist angle. 

These graphs can then be used as a preliminary orientation 
by the designer to achieve the optimal or near optimal 
combination of laminate configuration layer angle within the 
laminate that fit the specific requirements for the particular 
blade in focus. 

LITERATURE REVIEW 

This section reviews briefly some pertinent publications 
listed as References in the thesis report. Other related 
literature found during the scope of this research, is listed 
under Bibliography. The topics reviewed are classified under 
the following categories: 

. Finite Element Theory. 

. Composite Material Theory. 

. Vibrations. 

. Helicopter Theory. 

. Design Applications (Modeling). 

. Auxiliary Software. 

In the Finite Element Theory, relevant to present 
research, some basic references are: Batozf 19801 that 

describes the finite element procedures in engineering 
analysis using Discrete Kirchoff Theory ( DKT ) , Cookfl98l1 on 
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basic concepts and applications of finite element analysis, 
Craig[1981] describing computer methods in structural dynamics 
and f CASA/GIFTS, 1987] that is the reference manual for the 
finite element code used for the present (static and dynamic) 
work. Its graphical capabilities are very convenient to the 
analyst in a manner that simplifies the visualization of the 
effects of loads and couplings in the structure, making 
difficult concepts more accessible to understanding. 

Tsai and Pagano[1968] establish a notation in which the 
composite lamina properties are invariant with respect to the 
axis of rotation. This approach provides a very useful way to 
compare various material systems; the laminate theory is well 
documented, see for example Vinsonfl987] and Jonesfl975] both 
offer a good source on the mechanics and behavior of composite 
materials . 

Yntema[1955] developed a very useful tool to estimate 
bending frequencies of rotating beams in which the stiffness 
effect of RPM is shown for several blades (beam) 
configurations. Schi Ihansl f 1 958 ] uses another approach to 
estimate bending frequencies of rotating cantilever beams, 
introducing the effect of the angle made by the minor axis of 
inertia with the direction of the tangential circular 
velocity. Woodfl965] show^s a parametric investigation of 
aerodynamic and aeroelastic characteristics of articulated and 
hingeless rotor systems. Pr i tchard[ 1 988 ] presents an optimal 
placement of tuning masses to reduce vibration levels in 
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helicopter blades. Kottapall i [ 1 983 1 approaches the vibration 
reduction problem by modifying the blade torsional response. 
This is an innovative approach in that he modifies the 
aerodynamic response instead of adding mass, which is the 
usual conventional way to solve the vibration reduction 
problem . 

The helicopter theory in general is found comprehensively 
in Johnson [ 1 980 ] and is complemented in this research by 
Prouty[1986] in the understanding design trends. 

Hoskin[1986] is a general informative book in the 
aeronautical utilization of composite materials, and in 
Lake[1988] there is a preliminary investigation of finite 
element modeling of composite material rotor blades. 
Nixon[1987] addresses the extension-twist coupling of 
composite circular tubes applied to tilt rotor design. In 
McVeigh [1983] the aerodynamic design of XV-15 advanced 
composite tilt rotor blade requirements are reviev^ed. Finally 
Hodges[1987] presents a comparison of composite material rotor 
blade models, using two different methods of analysis: this 
model v;as chosen as a prime reference to build the basic model 
used throughout this work. 

Several software tools made it possible to manipulate the 
immense volume of data transforming, it into a suitable and 
compact presentation form, enabling the conclusions to be 
obtained. Kelly! 1988] was the basic reference to use Word 
Perfect 5.0 in the word processing and graphical editing jobs. 
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Holtf 19881 is a reference that helps to manage the data into 
Lotus 1-2-3 files and to get the input graph files to Word 
Perfect and to SURFER[ 1 987 1 . The SURFER program was used to 
obtain the 3-D carpet plots of extension-twist coupling 
analysis . 

C. THESIS OUTLINE 

This thesis report is divided into three areas, often 
found during a design process. 

It begins by modeling a rotor blade with laminate 
composite material. A numerical solution is obtained, 
exploiting the extension-twist coupling inherent when 
asymmetries are present in the laminate construction. 
Discrepancies encountered in the results compared with other 
rotor models are reported. These differences raised questions 
fostering the research towards the possibilities of different 
rotor constructions, starting with the same basic laminate. 

The different rotor configurations obtained are then 
subject to both static and dynamic analysis. The static 
analysis reveals the different responses of these models, 
while the dynamic behavior remains almost invariant, yielding 
only small variations in the first four fundamental vibration 
modes . 

Analysis of the results is done, and it is believed that 
a new insight into finite element analysis of composite 
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laminates is provided. This new insight is better noticed w^hen 



there are asymmetries present within the 
construction . 



laminate 
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THEORET I CAL EO RMU L AT I ON 



A • BASIC CONSIDERATIONS 

To understand how composite material works and its 
particular behavior that is useful in "designing" the 
material, a macroscopic approach is resorted to. 

Composite means that two or more materials are combined 
and usually preserve their best qualities, and, often exhibit 
properties that neither one possesses by itself. 

Some of the properties that can be improved by making a 
composite material are: 

. strength 
. stiffness 

. corrosion resistance 
. wear resistance 
. life in fatigue 
. thermal insulation 
. thermal conductivity 
. acoustical insulation 

It may be noted that not all these properties are improved 
at the same time, and usually it may not be required. 
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An Orthotropic^ material has properties which are 
different in three mutually perpendicular directions at a 
point, in other words has three different planes of symmetry. 

A lamina is the basic form of any laminated fiber- 
reinforced composite, it consists of a plane mounting of 
fibers imbedded in a matrix. 

In the macro-mechanical behavior of the lamina, only 
averaged apparent mechanical properties, or better, the 
stress-strain relations for an orthotropic material under 
plane stress conditions, transformed to directions not aligned 
with the principal directions of the lamina, are considered. 

This describes the behavior of the orthotropic material 
whose laminae have arbitrary directions, other than the 
natural geometric directions of the structure to be analyzed. 

B • STRES S — STRAIN RELATIONS 

The stress-strain relations or generalized Hooke’s law for 
anisotropic material will be reviewed briefly, and the basic 
restriction of linear elastic behavior assumption applies. 

The following derivations can be found in more detail in 
Vi ns on [ 1 987 ] . 

The elasticity tensor has nine independent elasticity 
constants, because of symmetry and the three mutually 



- Orthogonally anisotropic. 
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orthogonal planes of symmetry ( Ci 6 =C 26 =C 36 = C 45 = 0 ) ; 
that is, 
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and 



( 2 ) 



Oi = Cij X 



(i = l,2,3,4,5,6; j = l , 2 , 3 , 4 , 5 , 6 ) = 



In these relations the hygrothermal and thermal effects 
are not included for the sake of simplicity of this report. 

Through a force equilibrium study, analogous to the Mohr’s 
circle analysis [Jones, 1975] the stresses and strains in the 



^The material directions 1 to 6 are equivalent to 
cartesian directions: x , y , z , y z , xz , xy respectively in material 
axis. 
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material directions can be related to arbitrary^ directions 
using the following relations: 
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The transformation matrix is constructed by defining 
m=cos 6, n=sin 6; where 0 is the ply orientation angle, 
defined with respect to material and structural axes, assuming 
positive sign in counterclockwise direction (Figure 1). 



^ More 

geometry for 



convenient , 
example . 
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axis related 



to the structure 
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Figure 1 Positive rotation of principal material axes from 
arbitrary xy axes [ Jones , 1975 ] . 



The inverse of 
terms of direction 



transformation matrix 
cosines m and n as: 



[T] 



is expressed in 



(5) 
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For a lamina of composite material the modern notation 
uses Qjj for the quantities in the material stiffness matrix 
and are defined as: 
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If the lamina has same properties in both 2 and 3 
directions, then Vj 2 = v^ 3 , V 3 ^ = V 32 = 0 , Gi 2 “^i 3 ) £22 = ^33 and a 
simplification in equations ( 6 ) lead to equations (7). 
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To compare various material systems in the design of a 
composite structure, Tsai and Pagano[1969] arrived at some 
invariants with respect to the axis of rotation: 
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These invariants are representative of the lamina 
properties. The [Q] matrix may be transformed to x-y axes from 
1-2 axes by the usual coordinate transformation. 
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[Q] = [T]‘^ [Q][T] 



In terms of the invariants Us, the fQ] matrix elements 
are written as: 



U2Cos(26) + U2Cos(46) 



Q 22 = U2Cos( 26) + UgCos(46) 



( 10 



^' 4 “ U2Cos(46) 



Qcc= U.- U„cos(40) 
DO 5 3 



Qi 6= ^ U2sin(26) + U2sin(46) 



^26" 2 ~ U2sin(40) 
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After the definitions in (10), the stress-strain relations 
for the lamina of a N ply laminate (Figure 2) are given by: 



(11) 
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Where are defined earlier. 




Figure 2 Laminate construction, f 90°/+4 5°/0°/-45°/90° J 
characterized in this case to show the nomenclature for 
the lay-up. 
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LAMINATE ANALYS I S 



When two or more laminae are bonded together to obtain an 
integral structural component, the laminae principal material 
directions are established to produce a component that can 
resist loads in several directions (Figure 2). 

The following procedures enable to obtain the stiffness 
matrix of such a composite material from the lamina 
properties, using arbitrary combinations of materials and 
layer angles. 

For the lamina of a N layer laminate, equation (11) 
represents the stress-strain relations. 

In order to obtain the strain-displacement relations, the 
displacements are defined as follows: 





u(x,y,z) 


- u^(x,y) 


+ za(x, 


y) 


( 12 ) 


v(x,y,z) 


= v^(x,y) 


+ zf3(x, 


y) 




w(x, y , z ) 


- w(x,y) 







where u, v and w are displacements in x, y and z directions 
respectively 
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( 13 ) 



a 



and 



e = 



6x 



6y 



are the negative of the slope, related to the rotations. 



The mid-surface strains, then are given by 



( 14 ) 
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while the curvatures are given by 
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derivative . 



expressions 6 means partial 



Defining h as laminate thickness, (Figure 3) as the 

vectorial distance from the panel mid-plane, N, M, Q as stress 
resultants (Figure 4), stress couples and shear resultants 
respectively. Relations can be established relating these 
resultants to the strains and curvatures. This matrix is 
henceforth referred as the ABD matrix. 

This nomenclature is adopted for clarity and simplicity 
as may be noted from the expressions that follow. 

The stress resultants are expressed by integrating 
stresses across the thickness of the plate and is given by: 
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and the moment resultants take the following form 
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Figure 3 Nomenclature for the Stacking Sequence 
[Vinson , 1987 ] . 




Figure 4 Positive directions for Stress Resultants and 
Stress Couples for a Plate [Vinson, 1987] . 
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Relating resultants with the inplane strains and 
curvatures through expressions (11) through (14), matrix form 
is given by: 



[N] = [A].[6j + [B].[K1 

(18) 

[Ml = [B].[ej + [D].[K1 



In expression (18), the matrices A, B and D are called 
extensional stiffness matrix, bending-stretching coupling 
matrix and flexural stiffness matrix respectively, where: 
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Assuming that transverse shear stresses are distributed 
parabol ically (even though there exists discontinuities at 
the laminae interfaces), the shear resultants may be written 
as : 



( 20 ) 
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where 



(21 ) A. .= f , ( Q. . ). 

ij 4 k=l ij k 



^k ^k-1 3 * ^^k~ ^k-l'* ~2 

n 



i , j = 4,5 only 

The stresses and moments resultants can be written in a 
compact form that is very convenient to see the overall 
relations within a composite material laminate, specifically 
the ABD matrix. The ABD matrix is obtained as a combination 
of the expressions in (18). Expressions (18) and (22) are 
equivalent and in expression (22) the ABD matrix is defined 
explicitly. 
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ID . EFFECTS OF FLY ANGLE e 

To get an insight of the overall influence of the ply 
angle 6 on the laminate stiffness, it is useful to observe 
the effects of 6 for a given lamina. As each element of the 
ABD matrix, A^^, B^ j , and D^^, is a function of the material 
stiffness coefficient Qij, the layer thickness tj^ and mid-plane 
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vectorial distance the Qi^’s represent the stiffness 
variations due to ply angle changes in a given lamina. 

The property of the composite material that reveals its 
advantages over metallic or plastic is that it can be tailored 
to fit a given structural shape and a specific set of loads 
and operating conditions. 

The terms of the ABD matrix are functions of the material 
properties and the ply angle 6. If the properties are kept 
constant and within typical values of most fiber resin types 
of composite, it is possible to plot the values of versus 
0 and get a set of seven very useful graphs. These will permit 
an initial appreciation of such effects, and how they can be 
combined to meet a specific structural requirement. 

Figures 5 to 11 show how the elements of the "material 
stiffness" matrix ABD can be tailored by changing the ply 
angle, taking full advantage of this characteristic behavior 
of laminated composite material. 
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Figure 6 



22 



Qbar66 X THETA 




Figure 7 



Qbar12 < THETA 




Figure 8 
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Figure 9 




Figure 10 
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In figure 11 all Q’s are superimposed, on the same plot, 
to give a better perspective of the trends in doing this kind 
of layer design to fit specific structural constraints. When 
Ajg and A 26 are non-zero there is stretching-shearing coupling, 
while non-zero means twisting-stretching coupling, and non- 
zero B 2 s indicates bending-shearing coupling; bending- twist ing 
coupling is originated in non-zero values of and D 2 e • 

In normal design procedures, these terms are made zeroes 
by selecting appropriate stacking sequence; cross ply (0°/90° 
plies combination) or angle ply laminates {-0/+0 plies 
combination). Other important factor to be considered is the 
symmetry with respect to the mid-surface plane. 
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The stress strain relations or constitutive relations ( ABD 



matrix), combined with the proper equations of equilibrium and 
the strain-displacement relations form the basis for analyzing 
beams, plates or shells. 

These basic relations lead to the motivation and reasoning 
behind the model tailoring to fit a design requirement. 
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III. 



ROTOR, MODELING 



A. TYRICAL CURRENT ROTOR 

MODELING 



Current rotor design is a compromise between the different 
requirements in hover and forward flight, and the analytical 
models reflect that tendency or limitation. For instance, the 
optimum twist for' hover is not the same in forward flight for 
tilt rotors as well as for helicopters. Indeed, very often, 
the designers avoid extens ion-twist coupling w’ith laminated 
composite materials because of certain undesirable effects. 
The result is that one of the two flight modes will have less 
than optimal propulsive efficiency, depending on design phase 
decisions driven by requirements and/or aircraft mission 
prof iles . 

Vith the tilt rotor concept of mixing fixed wing and 
rotorcraft technology, these design trends were brought to 
surface once more. 

B . RROROSED SOLUTION 

Kixon[1987] proposed a solution for this sort of dilemma. 
The model used an extension-twist coupled rotor blade w'ith 
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100% RPM^ in hover and 80% RPM in forward flight. These are 
reasonable ratios between hover and forward flight for a tilt 
rotor type aircraft. The difference in RPM provides a change 
in centrifugal force, inducing a change in the twist angle 
distribution, which can be used to optimize the rotor for both 
regimes • 

The use of composite materials makes it possible to design 
a rotor blade that takes advantage of the changing loads, 
rendering the required twist distribution optimum in each 
flight mode. At this point one question arises: 

P/hat is the magnitude of twist available from a particular 

geometry and material system? 

The ansvv^er to that question lies in a method of 
prediction, with a model that behaves as close to the actual 
blade as possible yet inexpensive enough to permit a wide 
spectrum parametric study, providing a source to evaluate 
several geometry and material systems. Certainly such modeling 
should be verified by experimental data. 

This research is aimed to address the analytical part o f 
the above question and provide design tools to rotor blade 
systems . 



^ Reference value for 100% RPM is 217 rad/s. in that 
scale of rotor blade. 
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The torsional moments that affect a rotating blade can 
modify substantially the twist distribution. 

In this numeric model, the pitching moments generated by 
aerodynamic loading are not taken into account, the emphasis 
being on the inertia and extension-twist coupling as sources 
of torsional moments. 

An assumed lift distribution acts on the model, but 
compared with the inertia forces, its effects are negligible 
on the rotating blade. The inertia forces acting on the 
rotating blade produce extension which causes twist by 
structural coupling. However, the same inertia forces produce 
a centrifugal flattening [Prouty, 1986] effect called tennis 
racket effect (because of the tendency of a tennis racket to 
align its plane with the plane of rotation as it is swung in 
an arc ) . 

The forces acting in a helicopter blade can be seen in 
Figure 12, and it may be noted that in this case the effect 
is reduced with balancing weights (Chinese weights). 

C* GIFTS CAPABILITIES 

Gifts capabilities include generation of any type of 
structural model, with the following elements: 

. one-dimension rods or beams 
. two-dimension plates 
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. three-dimension shells, solid and complex stiffened 

shells . 

The elements can be selected from a library of options. 

The materials can be created by the user or selected from 
a library of standard definitions. The program allows the user 
to define isotropic materials and laminated composite 
materials, that can have up to a maximum of hundred layers. 

Loads can be concentrated or distributed, static as well 
as function of time. 

Static and dynamic analysis are performed on the models 
generated. The dynamic analysis provides free vibration 
response and forced response (transient computations), and 
model superposition . 

In using GIFTS as the finite element program to model the 
blade, the effects of the centrifugal force are taken into 
account with the forces acting on the lumped masses as shown 
in Figure 12, radially from the center of rotation. In using 
other finite element programs this feature must be verified 
for correct usage and application. The tennis racket effect 
can affect the twist angle to a large extent. If the 
centrifugal forces act in parallel lines instead of passing 
through the center of rotation, there will be no centrifugal 
flattening effect and results in larger twist angles. 

The basic rotor model selected to apply the proposed 
method is the one analyzed by Hodges f 1 98 7 1 , this model was 
analyzed using two different methods, namely by finite element 
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( MSC/NASTRAN ) and a 
coupled-beam model 

[Rehfield, 19851. The 
rotor geometry 
presented by 
Hodges[19871 was 
modeled in GIFTS using 
triangular TBS element 
( " flat triangular 
plate or shell 
element, including 

membrane, out of plane 
and fictitious in- 
plane bending 
stiffness . " [Batoz , 
19801). The results 
indicated a wide 
disparity of results. 
However, the present 
model appears to be 
within the assumption of linear and small displacements 
theory . 

The models in Figure 13 have a notation (M120 / Ml-20) 
that is described in the next section. A twist angle 
distribution discrepancy could be explained by the tennis 
racket effect probably not present in the model described in 




Figure 12 Twisting Moments Due to 
Centrifugal Forces f Prouty, 1986 ] . 
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Hodges f 1 987 ] , but the tip vertical displacement due to lift 
and weight with no rotation could not be explained by this 
argument. It may be observed that the displacements quoted in 
the referenced work are about 25% of the length, which appears 
to violate the linear theory that was used. 



BLADE DEFLECTION DUE TO LIFT AND WEIGHT 



COMPARtSON WITH MASA TM 0902^ 




BLAOe STATIOM c'»^CMeso 

O M120 MT- 20 c r> 4 ASA 



Figure 13 Comparison between GIFTS model and NASA TM 89024 
resul ts. 



This discrepancy 
models compared are 
effort was made to 



may be due to a possibility that the 
not exactly the same, although every 
take into account all the pertinent 
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geometric and material data reported. However this discrepancy 
was the driving factor that raised some questions, such as: 
What are the possibilities to model different structural 
configurations starting with the same number of layers, 
same lay-up design, same thickness and material 
properties? (perhaps by inadvertent oversight by the 
designer, analyst using canned programs, or manufacturing 
engineer . ) 

The answer to this was based on the following assumptions: 

. The structure is a D shape spar. 

. It is built using three shells, top, bottom and rear 
shell surfaces . 

. The top and bottom meet at the leading edge and the rear 
siiell closes the box. 

. Once established, the internal lay-up in the stacking 
sequence i-emains constant. 

With tliose assumptions a laminate to constitute one shell 
suid’ace may be destined, using only the Laminate Pi-incipal 
Orientation ( LPO) angle. For example, in the following 
laminate [ + 20°/ - 70°/ + 2 0°/ - 7 0^/ - 7 0°/ + 20° ] , the top layer angle 
is denoted +20° and is a six layer asymmetric laminate. 

Rotating the whole laminate 40° clockwise gives 
[ - 2 0° / + 70°/ - 2 0° / + 7 0° / + 7 0°/- 2 0° ] , whose Lami nate Principal 

Orientation angle is denoted as -20°. 

This laminate has i n t en t i onal 1 \ designed asymmetry tu gel 
the f^x t ens j on- 1 w i s t couiiing. Associated with this laminate 

O v-> 



is a normal vector whose orientation is seen in Figure 14, The 
D-spar then is built using three such laminates, with each 
laminate mounted with the normal vector facing inward or 
outward. This method of construction results in six different 
configurations of the rotor. These configurations and other 
variations within each one will be developed in the next 
section . 

D • VARIOUS ROTOR DESIGNS OR 
CONFIGURATIONS 

The six different configurations are designated Ml to M6 
and the LPO^ , varying from -90® to +90® in steps of 20®, 
provides the suffix to complete the model designation Ml-90^, 
M340, M6-60 etc. In a total of 11 LPO in each model. 

Different configurations may be visualized with the help 
of the scheme shown in Figure 14. 

Model I and IV can be built with only one laminate vrrapped 
around, one inverse of the other. 

In the program GIFTS, the outward vector for each surface 
is defined during the "GRID” definition, e.g. in a four sided 
grid, each shell is defined within a grid. The local X axis 



^ Laminate Principal Orientation Angle, referred to the 
top layer of the laminate used to construct that particular 
model . 

^ Ml-90 meaning model 1 with LPO = -90®; M340 meaning 
model 3 with LPO = 40®. 
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is oriented parallel to the first line of the grid, the local 
Z axis is the cross product of the first and second "line” in 
the grid definition, while the local Y axis is the cross 
product of the Z and X axes, as is usual in cartesian systems. 

The outward vector is oriented in a direction defined by 
the local Z axis. Layer angles are positive in the counter- 
clockwise direction in the XY plane for shell elements, as is 
the convention for laminated composite materials. 



+20 


a 


-70 




+20 




-70 




-70 




+20 


) 



END VIEW OF D SPAR ( 3 SHELLS) 






MODEL^ 1 



MODEL 2 




MODEL H 




NORMAL VECTOR ASSOCIATED TO 
LAMINATE STACKING SEQUENCE 



V 

MODEL 3 



MODEL 5 MODEL' 6 



Figure 14 Normal to surface vector associated with the 
layer sequence of the basic laminate and its position in 
relation to end view of the blade. 



It may be observed that, when an asymmetry is present in 
the basic laminate (material definition), the structural 



35 



response is highly dependent on how the grid and consequently 
the lay-up is defined. In the case of a symmetric laminate or 
isotropic material, this fact is never brought to surface 
because of its irrelevancy. However, dealing with asymmetric 
laminated composite materials, the analyst must always have 
this in mind. It may be noted that, in most cases the laminate 
is chosen to be symmetric. 

Table 1 gives in a compact form the description of various 
laminate constructions described. Included in that table are 
the line names actually employed in the input file M$**.SRC 
[Annex 1] {$ = 1,2, 3, 4, 5, 6 models; = top layer angle as 
defined previously, in the range -90 to +90), to develop the 
models generated in GIFTS. By modifying two areas in this 
input file, six models are generated that are considered under 
the present study, this also provides variation of the layer 
angle witliin each model. The actual modifications employed in 
the input file are given in LAYUP. AUX [Annex 2]. For each 
model (one of the six combinations obtainable by the change 
in the laminate associated vector, as given in Figure 14), 
the whole laminate was rotated in increments of 10°, 
referenced to the top layer angle (LPO). For each of these 
models^ the RPM , and consequently the centrifugal force, was 
varied in the range (130 rad/s to 300 rad/s). 

The model geometric characteristics are the same as in 
Hodges [ 1987 1 , a graphite epoxy composite rotor blade 
constructed out of Hercules IM6 fiber with Ciba-Geigy R6376 
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resin. The cured ply thickness is assumed to be 0.0055 in. The 
D-spar has a 35.23 in. radius, and begins at station 5.32 in., 
which has all degrees of freedom suppressed to emulate a rigid 
rotor connection to the hub. The rotor has a constant cross- 
section, defined by 13 nodes and is divided into 21 span-wise 
(3.0 in apart) stations. Models with finer meshes were used 
in the early stages to check for convergence. The mesh 
selected seemed to be a reasonable compromise between the 
computational time and accuracy. 

The orthotropic material properties of the lamina are 
listed below; 

Ejj , psi 23.1 X 10® 

E 22 , psi 1.4 X 10® 

Vj2 0.338 

Gi2,psi 0.73 X 10® 
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TABLE 1 - MODEL DEFINITION: 



Example 
Vector 
+ 20 k 
-70 
+ 20 
-70 
-70 
+ 20 



GRIDS, POINTS AND LINES 
5 L45 4 L14 1 L51 5 





G3 


G1 


G2 


L512 


L411 


L18 


L512 



12 L1112 11 L811 8 L128 12 

(not to scale, e.g. line L512 is 30.0 in length, 
L45 is ,286 in, and L14 is about 1.15 in) 

♦ = vector inward ; o = vector outward. 



POSSIBLE RIGHT HAND RULE NORMAL VECTOR TO SURFACE 

MODELS GRID DEFINITION 



Ml$$ 


VECTOR 


* 


* 


* 


GRID 

DEF. 


L512,L1112, 
L411 , L45// 


L411 ,L811 ,L18, 
L14// 


L18,L128,L512, 

L51// 


M2$$ 


VECTOR 


* 


* 


o 


GRID 

DEF. 


L512 , L1112 , 
L411 , L45// 


L411 ,L811 ,L18, 
L14// 


L512 ,L128 ,L18, 
L51// 


M3$$ 


VECTOR 


* 


o 


o 


GRID 

DEF. 


L512,L1112, 
L411 , L45// 


L18 , L811 , L411 , 
L14// 


L512,L128,L18, 

L51// 


M4$$ 


VECTOR 


o 


o 


o 


GRID 

DEF. 


L411 ,L1112, 
L512,L45// 


L18,L811 ,L411 , 
L14// 


L512 , L128 , L18 , 
L51// 


M5$$ 


VECTOR 


o 


o 


* 


GRID 

DEF. 


L411 ,L1112, 
L512,L45// 


L18,L811 ,L411 , 
L14// 


L18,L128.L512, 

L51// 


M6$$ 


VECTOR 


o 


* 


* 


GRID 

DEF. 


L411 ,L1112, 
L512,L45// 


L411 ,L811 ,L18, 
L14// 


L18,L128,L512, 

L51// 
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IV 



STATIC AND DYNAMIC 



ANALYS I S 

The six models described in chapter III are basically 

subjected to three main types of independent static loads. 
Further, a modal analysis is also performed to obtain the 
first four fundamental modes of vibration frequencies and mode 
shapes . 

The first load applied is a torque at the tip, not 

combined with any other load. 

The second load case consists of combined lift and blade 
weight at zero RPM , i.e. with no centrifugal forces. 

The third load case is the centrifugal force due to 

rotation of the blade. The lift and blade weight are present 
and the RPM varies from 130 rad/s to 300 rad/s. 

A. TORSIONAL LOAD 

A moment of 26.4 in-lb is applied distributed among the 
nodes at the tip as concentrated moments. 

This load case reveals the influence of the laminate 

principal orientation angle ( LPO ) , and the torsional stiffness 
for each model. The different response of the six different 
models also is presented for comparison purposes. 
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The results are presented in Figure 15 for the six models 
for LPO’s of +20° and -20° in each model. 

As much as 40% variation in the response may be observed 
in the comparisons. 



ANGULAR DISPLACEMENT COMPARISON 
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Figure 15 Angular Displacement at the Tip Due to Pure 
Torque. 



B. LIFT AND BLADE WEIGHT LOADS 

This loading is designed to get the combined effect of the 
lift and blade weight given by Hodgesf 1987 ] , expressed in 
terms of radial distance of the center as follows: 

(17) 1 = (0.02222 Ib/in) x r - 0.0123 Ib/in 
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This load is applied vertically in two spanwise rows at 
the leading and trailing edges, in a proportion that produces 
zero moments at the quarter chord. This load distribution is 
more homogeneous and induces less cross-section distortion, 
as opposed to the application along a line on upper surface 
as in beam type elements. This distributed line loading 
produces small torsional displacement ( lO*** degree). 

The results are presented in Figure 16 for the six models, 
with 2 values of LPO (+20® and -20®) for each model. 




Figure 16 Displacement at the Tip Due to Vertical Loads 



Along the Blade. 
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CENT Pi I FUGAL LOAD 



Next loading case investigated reveals the effect of 
varying centrifugal load, combined with changing the LPO angle 
within each model. The response yields twist angle induced by 
extension- twist coupling present in the structure due to 
asymmetric stacking sequence. 

The LPO, referred to the top layer angle, varies in the 
range of +/- 90°, the internal lay-up in the stacking sequence 
remaining constant, as previously defined. 

After some experiments to localize inflection points in 
the twist angle variation, the LPO to be investigated were 
chosen in the range of +/- 90°, giving a total of 11 laminate 
principal orientation angle for each of the six models. 




Figure 17 Overall Range of the Twist Angle within the Six 
Mode 1 s . 
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The loads due to lift and blade weight are present but 



their effects are observed to be negligible compared w 
those due to centrifugal loads. 

The RPM is varied from 130 rad/s to 300 rad/s 
increments of 10 rad/s. in a total of 19 load cases, 
combination of load cases and LPO give a total of 201 poi 
of twist angle for each of the six models. These points 
presented in six 3-D graphs to better appreciate the choi 
offered in this design approach. Those graphs present 
trends of the twist angle with RPM and LPO variations, 
twist angle range for each model is in Figure 17, while 
overall parametric information are presented in Figures 
throuoli 23. 
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Figure 18 Model One, Twist Angle Induced by Rotational 
Inertia Loads Combined With Laminate Principal Orientation 
( LPO ) , Referred to Top Layer Angle, 
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Figure 19 Model Two, Twist Angle Induced by Rotational 
Inertia Loads Combined With Laminate Principal Orientation 
( LPO ) , Referred to Top Layer Angle. 




Figure 20 Model Three, Twist Angle Induced by Rotational 
Inertia Loads Combined With Laminate Principal Orientation 
( LPO ) , Referred to Top Layer Angle. 
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Figure 21 Model Four, Twist Angle Induced by Rotational 
Inertia Loads Combined With Laminate Principal Orientation 
( LPO ) , Referred to Top Layer Angle. 
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MODEL 5 




Figure 22 Model Five, Twist Angle Induced by Rotational 
Inertia Loads Combined With Laminate Principal Orientation 
( LPO ) , Referred to Top Layer Angle. 
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Figure 23 Model Six, 
Inertia Loads Combined 
( LPO ) , Referred to Top 



Twist Angle Induced by Rotational 
With Laminate Principal Orientation 
Layer Angle. 
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MODAL ANALYSIS 



D . 

To observe the effects of each configuration (Model) and 
laminate principal orientation on the first four fundamental 
frequencies, a free vibration analysis is performed. 

The results of the modal® analysis are depicted in Figures 
24 and 25 for the six models and for LPO’s of -20° and +20° 
respectively . 




Figure 24 Frequency Variation among Different Models. 



® 1 FW means First Flapwise; 1 CW means First Chordwise 
and so on. 
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It may be noted that, even for simple models, this 
solution procedure is a time consuming task, only recommended 
for final or quasi final configurations in the design process, 
or using a very crude mesh lay-out. 

It can be seen that, for each vibration mode, the 
frequency remains practically constant. This reveals a low 
sensitivity of the frequency to changes in model 
configurations, considering the first four fundamental modes 
and same LPO’s. 




Figure 25 



Frequency 



Variation 



among Different 



Models. 
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V 



ANALYSIS OF RESULTS 



In this chapter the results obtained for each model will 
be presented and several observations are inferred. 

The present research was conducted to get an insight of 
the trends implicated into the two questions posed in Chapter 
III, sections B and C. This study also provides some insight 
into answering the questions posed in sections B and C of 
Chapter III. 



What is the magnitude of twist available from a particular 
geometry and material system? 



What are the possibilities to model different structural 
configurations starting with the same number of layers, 
same lay-up design, same thickness and material 
properties? 



The first question is a design approach to structural 
analysis. That approach will be elaborated in this chapter. 
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The second question opened a vast field for research 
activity that will be addressed in next Chapter. 

It can be seen in Figure 14 that the models 1 through 6 
are obtained simply by changing one of the normal vectors 
associated with the laminate, (in one of the three faces of 
the spar). The effects of these changes in the structural 
behavior of the blade may be divided into two main categories, 
static and dynamic behavior. 

A. DYNAMIC BEHAVIOR 

The free vibration characteristics are shown in Figures 
24 and 25. Figure 24 shows the frequencies for LPO of -20° and 
Figure 25 presents for LPO of 20°. The first four fundamental 
modes of vibration have very little sensitivity to the model 
type, and may be observed that the frequencies do not 
change, for the same or symmetric LPO (positive/negative). 

This behavior can be considered benign, or in other words 
auspicious, in a connotation that, once the designer arrived 
at the desired frequency values for the RPM range considered, 
the LPO or the model can be changed until a desired twist- 
extension coupling is achieved. 

These iterations can be done without looking very close 
to resonance problems, which must be checked anyway when the 



final model is to be selected. 



This assumption is made based on the validity of 
Yntema[1955] approach for these models. It is possible that 
the change in the stiffness matrix due to RPM affects the 
rotor built of laminated composite materials in a different 
way than isotropic rotors. Further research needs to be done 
along this direction to obtain comparable results for 
composite rotors. 

In Figures 24 or 25, each column corresponds to one mode 
of vibration at zero RPM and is a point in the 3-D plot for 
that specific mode (flatwise, chordw’ise or pitchwise^). The 
RPM and LPO may be v^aried as in the static case yielding 
carpet plots. The resulting plots would validate or not the 
assumption of Yntema*s values for composite blades. 

B- STATIC BEHAVIOR 

The static behavior is presented in Figures 18 through 23. 
The parameters considered are RPM and LPO as inputs and the 
twist-extension coupling manifested through the twist angle 
at the tip, as the output. This angle will be referred as TTA 
(tip twist angle) in the subsequent discussion. 

The TTA can be used as a measure of the twist angle ratio 
along the blade span. The twist angle is not a linear function 
of the radius but may be represented by a second degree 

^ Pitchvv’ise mode is a higher mode in this case and is not 
considered in this report. 
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polynomial, such that the TTA is a good parametric factor to 
infer the twist angle magnitude along the blade. 

As mentioned before, the six models are obtained by 
changing the orientation of the normal vector associated to 
the laminate, in the three shells that constitute the spar, 
one at a time as shown in Figure 14. 

In that figure it can be verified that Model 1 and Model 
4 can be built with only one unbroken piece of laminate, 
wrapped around a mold. Further, one model is just the inverse 
of the other, in other words Model 1 is the inside out of 
Model 4. 

The models 2 and 5 can also be considered as some sort of 
inverse, in a broader* sense; similar reasoning may be extended 
to models 3 and 6. 

With the abo\*e considerations, the six models can be 
analyzed in pairs, which is vei*y convenient from the point of 
view of getting the information in a simple and concise form, 
so tlial the designer* can quickly evaluate the model, keeping 
in \ iew the design requirements. 

The second important consequence of this inverse concept 
is that, Model 1 and Model 4 cari supply a spectrum of twist 
angle requirements, for a given loading condition. 

The importance of this reside in the manufacturing 
pr*ocess. in tliat both models may be built by filament wound 
techniqie-. Thi:’ process 



is less expensive and irore i*e liable 



as it eliminates the filament discontinuity in the shell’s 
junction, and requires less accurate quality control. 

The convention used in the modeling process has as a 
aftermath that positive angles have active aerodynamic 
stabilizing effect, acting as a propeller governor device. 

1 . Models 1 and 4 

Figures 18 and 24 exhibit the behavior of these two 

model s . 

The models, which were referred as inverse to each other, 
also show the response characteristics and behave as expected. 

Tliis may be noted by observing the displacement for a 
given RPM and Layer Principal Orientation Angle in Figure 26. 

2. Models 2 and 5 

Those models have the same inverse relationship in 
the TTA value. The main difference is that these two models 
are stiffei' than the other four, as shown in Figure 17. Figure 
27 shows these two models drawn together. 

3 . Models 3 and 6 

The inverse relationship also holds here and the 
stillness values are higher than that found in models 1 and 
A and lower than models 2 and 5 (Figure 17). See Figure 28 



loi this comparison. 



MODEL 1 








Figure 26 Models 1 and 4 static behavior comparison, 
are RPM and Laminate Principal Orientation ( LPO ) 
output is Tip Twist Angle ( TTA ) . 



inputs 
angl e . 
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MODEL 2 







^ U ^ 



MODEL 5 



>x\ 




Figure 27 Models 2 and 5 static behavior comparison, inputs 
are RPM and Laminate Principal Orientation ( LPO ) angle, 
output is Tip Twist Angle ( TTA ) , 
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f-lODEL 3 




Figure 28 Models 3 and 6 static behavior comparison , inputs 
are RPM and Laminate Principal Orientation ( LPO ) angle, 
output is Tip Twist Angle ( TTA } . 
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VI 



CONCLUSIONS AND SCOI>E 



FOR FUTURE RESEARCH 
A • CONCLUSIONS 

This report presents a detailed analysis of a typical 
rotor blade, the focus being the static behavior in the 
presence of extension-twist coupling due to asymmetric 
stacking sequence of laminas. 

Different extension-twist coupling effects are obtained, 
by changing the ply orientation with respect to the body axes, 
while keeping the stacking sequence, thickness and internal 
lay-up sequence as invariants. 

These effects are presented in terms of carpet plots. The 
extension-twist coupling is measured through the angular 
displacement at the tip of the blade as function of RPM and 
Laminate Principal Orientation ( LPO ) angle variations. 

Free vibration analysis is performed for two LPO ’ s of the 
six different models. The first four fundamental frequencies 
show small sensitivity to these specific structural changes. 

The static analysis revealed a sort of ”ant isymmetr ic” 
behavior within the models, allowing them to be grouped in 
pairs. Such behavior can be of use in the design process, to 
trim structural response to prescribed loads. 
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A structural analyst should be very careful in dealing 
with composite laminated structures, especially if the analyst 
is experienced only in structural analysis of isotropic 
materials. It was noted that it is very easy to be mistaken, 
mainly if there are asymmetries present within the laminate. 
It is very important to know how the local coordinates systems 
are established in the finite element program being used. 
These axes will determine the relative position of the 
laminate with respect to the structure. If asymmetries are 
present, it is possible to get different structural response 
and have different structural configurations, all starting 
with the same basic laminate. 

B • SCOBE EOE EUTUEE RESEARCH 



Future research may be beneficial into two main areas. 

The first one being to validate these models with 
experimental data, providing the necessary confidence in 
further utilization of the method envisaged. 

Once the validation is done, the method can be corrected 
and trimmed for operational use. 

Then a field that calls for attention is a dynamic 
characteristics mapping, that is, with the RPM and LPO as 
parameters obtain the fundamental frequencies as output. 

This would be a 3-D composite material extension of 
Yntema’s [1955] work for beams of isotropic material. 
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ANNEX 1 



M 1 2 O . S ¥tC 



t DEFINES BLADE MODEL WITH NO ELIPSE PROJECTION 

$ ( LOST 3.976 DIGITS IN DECOM ) 

S- ( ~ 20 MIN TO RUN 2 LDCASES ) 

$ NO ROOT DEFINED JUST SUPRESSED DOF OF FIRST RIB 

$ 

KPOINT 

1/ .65,5.23,0 
2/ .455,5.23, . 1 103 
3/0,5.23, . 1545 
4/- .49,5.23, . 1433 
5/-. 49, 5. 23, - . 1433 
6/0,5.23, - . 1545 
7/ . 455,5.23, -.1103 
8/. 65, 35. 23,0 
9/ .455,35.23, .1103 
10/0,35.23, . 1545 
11/-. 49, 35. 23, . 1433 
12/- . 49,35. 23, - . 1433 
13/0,35.23,-. 15^5 
14/ . 455, 35 , 23, - . 1 103// 

S 

$ LINES GENERATION 

S 

SL INE 
L34/3,4,2 
L45/4,5,4 
L56/5 .6,2 
L18/1 ,8,21 
L4 1 1/4 , 1 1 , 21 
L512/5, 12,21 
LlOll/10,11,2 
LI 1 12/1 1 , 12,4 
L1213/12, 13,2// 

$ 

$ ARC GENERATION 

$ 

CARC 

Cl 3/1 ,2,3,5 
C61 /6 , 7 , 1 ,5 
CSlO/8,9. 10,5 
C138/13, 14,8,5// 
f 

4 COMPOSITE LINES GENERATION 

5 

COMP;. INE 

L14/C13.L3- 

L5 1 / '_ 5o , C6 I 

lBI 1 '381 _ 101 1 

LI 26/^ 1213,3138// 



$ BLADE MATERIAL AND ELEMENTS DEFINITION (IM6/R6376) 

$ BLADE MASS DISTRIBUTION = .0123 LB/IN = .00038233 SLUG/IN 

S (COMPOSITE DENSITY: .001796 SLUG/IN"3) 

$ (TUNGSTEN DENSITY =.021116 SLUG/IN"3) 

$ (MEAN DENSITY =.004418 SLUG/IN*3) 

$ (COMPOSITE AREA =.086525 IN*2) 

$ (COMPOSITE VOLUME =2.5957 IN‘3) 

$ (TUNGSTEN AREA =.01358 IN*2 => DIA=.1315 IN) 

# (MODEL FOCUSED LENGHT=30 IN; TOTAL LENGHT=35.23 IN) 

ORMAT ,53 

1/5.07E5,23. 1E6, .338, .73E6, .004418 
8. 12E3, 1 .4E6, 1 .41E4// 

$ 

$ LAYER DEFINITION [+20/ -70/+20/- 70/ - 70/+20 ] 

$ THICKNESS = .0055 IN/LAYER * 6 = .033 IN 

COMPTH 

1 , 6/1 , 1 , 1 , 1 , 1 , 1 
.0055, .01375,20 
.0055, .00825,-70 
.0055, .00275,20 
.0055, - .00275, -70 
.0055, - .00825, -70 
.0055, - .01375,20// 

S 

$ ELEMENTS 

$ GRIDS MUST BE GENERATED IN SAME ORIENTATION 

S STABLISH A VIEW POINT NORMAL TO SURFACE OUTSIDE (OR INSIDE) 

X WITHOUT CHANGING THIS REFERENCE GENERATE THE GRIDS FOR 

$ ORTHOTROPIC MATERIALS (FILAMENT WOUNDED THINKING) 

$ 

f THE FIRST SIDE OF GRID DETERMINES THE REFERENCE FOR MATANG 

$ 

GETY/TB3/ 1,1// 

GRID4 

G1/L41 1 , LSI 1 , L 18 , L 14 
G2/L18,L128,L512,L51 
G3/L512, L 1 1 12, L41 1 , L45// 

END 



AlsTNEX: 2 



LAYUE> . AUX 



FILE TO BE USED IN BUILDING THE SEVERAL LAYUPS IN 
EACH M$**.SRC FILES 

f r layup number 

** = ANGLE OF TOP LAYER WITH RESPECT TO NORMAL VETOR DIRECTION 

1- LAYUP GRID DEFINITION -- 6 MODELS -- 

M1**.SRC 

GRID4 

G1/L41 1 , LSI 1 , L 18, L 14 
G2/Lie,L128,L5l2,L51 
G3/L512,L1112,L411 ,L45// 



M2*» .SRC 
GRID4 

G1/L41 1 ,L811 ,L18,L14 
G2/L512,L128,L18,L51 
G3/L512 , L 1 1 12 , L4 1 1 , L45// 

M3** . SRC 
GRID4 

Gl/Lie,L811,L411,L14 

G2/L512, Lies, Lie, L51 
G3/L512 , L 1 1 12 , L4 1 1 , L45// 

M4** . SRC 
GRID^ 

Gl/LIB, LSI 1 . L^l 1 ,L14 
G2/L512 , L 128 , L 18, LSI 
G3/L^l 1 , L 1 1 12, LS12, L^S// 

MS** . SRC 
GRID4 

Gl/L 18 , LSI 1 , L4 1 1 , L 14 
G2/L18,^128,L512,LS1 
G3/L41 1 , L 1 1 12 , LSI 2 , L45// 

M6** .SRC 
GR : D^ 

G1/L411,lS 11,L1S,L14 
G2/LlS,_l2S,LSl2,L51 
G3/L41 1 , . 1 1 12 , LS12 , .4S// 



2- TOP LAYER ANGLE -- 18 ANGLES -- 

M$00. SRC 
.0055, .01375,0 
.0055, .00825,-90 
.0055, .00275,0 
.0055, -.00275, -90 
.0055, - .00825, -90 
.0055, -.01575,0// 

♦ ♦*************-)Mt-**-X-***************->M(--tMMMM»-** 

MSIO. SRC 
. 0055 , .01375, 10 
.0055, .00825,-80 
.0055, .00275, 10 
. 0055 , - . 00275 , -80 
. 0055 , - . 00825 , -80 
.0055, - .01375, 10// 

N$20.5RC 
.0055, .01375,20 
. 0055 , . 00825 , -70 
. 0055 , . 00Z75 , 20 
.0055, - .002~5, -70 
.0055,-.00825,-70 
.0055, -.01 375. 20/'/ 

NJ30. 5RC 
. 0055 , .01 375 , 30 
. 0055 , . 00825 . - 60 
.0055, .00275,30 
. 0055 , - . 00275 , - 60 
. 0055 , - . 00825 , -60 
.0055, - .01375,30// 

NSAO. SRC 
.0055, .01375,40 
.0055, .00825,-50 
.0055, .00275,40 
.0055, - .00275, -50 
.0055, - .00825, -50 
. 0055 .-.01 375 , 40/ / 



MJ50.SRC 
. 0055 , .01375,50 
.0055, .00825, -40 
.0055, .00275,50 
.0055, -.00275, -40 
.0055, - .00825,-40 
.0055,- .01375,50// 

M$60.SRC 
.0055, .01375,60 
.0055, .00825,-30 
.0055, .00275,60 
.0055, - .00275, -30 
.0055. - .00825, -30 
.0055, - .01375,60// 

M$70.SRC 
.0055, .01375,70 
.0055, .00825,-20 
.0055, .00275,70 
.0055, - .00275,-20 
. 0055 , - . 00825 , - 20 
.0055, - .01375,70// 

Mi'SC' . 5RC 
. 0055 , . 0 1 375 . 80 
. 0055 , . 00625 ,-10 
.0055, .00275,80 
.0055. - .00275 , - 10 
.0055, - .00825,-10 
.0055, - .01375.90// 



90 . 5RC 



0055, 


.01375,90 


0055, 


. 00625 , 0 


0055 , 


.00275,90 


0055, 


- .00275,0 


0055 . 


- . 00825 , 0 


0055 . 


- .013~5.9 



M$-90. 5RZ 



COD':-, 


. 01375 .-'50 


0055 , 


. 00825, C' 


0:’55 , 


.ooz:’5, - = ■ 


'*/ ^ “ 


- . o: 275 . 0 


OC ”5 » 


- . I C'625 , 


1^' , cr 


_ ~ “7 ~ _ q 



“ BO . SRL 

.0055, .01375,-80 
.0055, .00825. 10 
.0055, .00275,-80 
.0055, - .00275, 10 
.0055, -.00825, 10 
.0055, - .01375, -80// 

MS-70. SRC 
.0055, .01375,-70 
.0055, .00825,20 
.0055, .00275,-70 
.0055, - .00275,20 
.0055, - .00825,20 
.0055, - .01375, -70// 

MS-60. SRC 
.0055, .01375,-60 
.0055, .00825,30 
. 0055 , . 00275 , -60 
.0055, -.00275, 30 
.0055. - .00825,20 
. 0055 ,-.01 375 , -60// 

♦ ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦•♦-ini. 



MS -5'/ . SRC 



0055, 


. 0 1 w 7 5 , 


-5C' 


0055, 


.00325, 


40 


0055, 


. 00275 , 


-50 


0055 , • 


- . 00275 




0055, • 


- .00825 


,^C 


0055, ■ 


- .01375 


_ cr 

1 



MS -40 . SRC 
.0055, .01375,-40 
.0055, .00825,50 
.0055, .00275,-40 
.0055, - .00275,50 
.0055, - .00825,50 
.0055, -.01 375, -40// 



MS -30. 


SRC 


. 0055 , 


.01375 , -30 


.0055. 


. 00825 , 60 


. 0055 , 


. 00275 , -30 


. 0055 . 


- .0!0275.t0 


. 0055 , 


- . 00825 , tC 


. 0055 . 


- . 0:375 , -3 






M$-20.SRC 
.0055, .01575,-20 
.0055, .00825,70 
.0055, .00275,-20 
.0055, - .00275,70 
.0055, - .00825,70 
.0055, - .01375,-20// 

MS-10. SRC 



.0055, .01375,-10 
.0055, .00825,80 
.0055, .00275,-10 
. 0055 , - . 00275 , 80 
.0055, - .00325,80 
.0055, -.01375, -10// 
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